Nature of the metal-insulator transition in the half— filled t — t' Hubbard chain 
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We study the quantum phase transition from an insulator to a metal realized at t' = t'^ > 
0.5t in the ground state of the half-filled Hubbard chain with both nearest-neighbor (t) and next- 
nearest-neighbor {t') hopping. The study is carried out using the bosonization approach and density 
matrix renormalization group calculations. An effective low-energy Hamiltonian that describes the 
insulator-metal transition is derived. We find that the gross features of the phase diagram are well- 
described by the standard theory of commensurate-incommensurate transitions in a wide range of 
parameters. We also obtain an analytical expression for the insulator-metal transition line t'i.{U,t). 
We argue that close to the insulator-metal transition line, a crossover to a regime corresponding to 
an infinite-order transition takes place. We present results of density-matrix-renormalization-group 
calculations of spin and charge distribution in various sectors of the phase diagram. The numerical 
results support the picture derived from the efi'ective theory and give evidence for the complete 
separation of the transitions involving spin and charge degrees of freedom. 

PACS numbers: 71.27.+a, 71.30.-fh 



I. INTRODUCTION 

During the last decades, the Mott metal-insulator tran- 
sition has been the subject of great interesti ^i^i^ In the 
canonical model for this transition - the single-band Hub- 
bard model - the origin of the insulating behavior is 
the on-site Coulomb repulsion between electrons. For 
an average density of one electron per site, the transition 
from the metallic to the insulating phase is expected to 
occur when the electron-electron interaction strength U 
is of the order of the delocalization energy (which is a 
few times the hopping amplitude t). The critical value 
(U /t)c turns out to be quite independent of the specific 
band structurei^ It is important to recall that the Mott 
transition is often preceded by antiferromagnetic order- 
ing, which usually leads to insulating behavior and thus 
masks the Mott phenomenon. 

While the underlying mechanism driving the Mott 
transition is by now well understood, many questions re- 
main open, especially about the region close to the tran- 
sition point where perturbative approaches fail to pro- 
vide reliable answers. The situation is more fortunate in 
one dimension, where non-perturbative analytical meth- 
ods together with well-controlled numerical approaches 
allow in many cases to determine both the ground state 
and the low-lying excited statesi^*Sii However, even in 
one dimension, apart from the exactly solvable cases, a 
full treatment of the fundamental issues related to the 
Mott transition still constitutes a hard and long-standing 
problem. 

In this paper, we study the t — t' Hubbard chain which 
includes both nearest (t) and next-nearest-neighbor (t') 
hopping terms. We limit ourselves to an average density 
of one electron per site (the half-filled band case). De- 
pending on the ratio between t' and t, the system has two 
or four Fermi points. Correspondingly, it shows one- or 



two-band behavior and has a rich phase diagram. There- 
fore, it is not surprising that the model has been the 
subject of intensive analytical and numerical studies, in- 
cluding a weak-coupling renormalization group analysis;^ 
DMRG calculations for charge and spin gaps^iiSiiiiiSii^, 
the electric susceptibilityji^ the momentum distribution 
function^iSiiSiii and the conductivitjii^ as well as, very 
recently, a variational technique^ 

Unfortunately, conflicting results have been reported 
for the transition region, in particular regarding the char- 
acter of the transition, the number of different phases and 
the number of gapless modes. In this paper, we hope to 
settle some of the unresolved issues using a combined 
analytical and numerical analysis. We focus our atten- 
tion on the metal-insulator transition as a function of t' 
for a fixed on-site repulsion U. An effective continuum 
theory allows us to show that in the parameter range 
0.5t < t' < t'^ the system exhibits the characteristic be- 
havior of a commensurate-incommensurate transition^ 
Close to the transition point, additional scattering pro- 
cesses characteristic of two-band behavior set iuiS We ar- 
gue that these processes induce a crossover to Kosterlitz- 
Thouless type critical behavior, as found in Ref. Q. Ad- 
ditional support for this conclusion is provided by DMRG 
calculations for chains of up to L = 128 sites and large 
U . The numerical analysis also allows us to study the 
gaps in the excitation spectrum as well as the charge and 
spin density distributions. 

The paper is organized as follows. In Section II, gen- 
eral properties of the model including the strong coupling 
limit and the structure of the phase diagram are briefiy 
reviewed. Section HI shows that the weak-coupling 
bosonization approach leads to the quantum sine-Gordon 
field theory, the standard model for commensurate- 
incommensurate transitions. In Section IV, the results 
of the numerical analysis are presented, including the ex- 
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citation spectrum and density distributions. Special at- 
tention is given to the regime of large t' where the model 
represents a system of two coupled chains. The results 
are summarized in Section V. 



II. THE t t' HUBBARD CHAIN 

The one-dimcnsional t — t' Hubbard model is defined 
by the Hamiltonian 

+ [/ ^ {n,.^ - 1/2) K- X - 1/2) , (1) 



where c j ^ ('^jo-) ^"^^ electron creation (annihilation) op- 
erators on site j and, with spin projection a =t,i, 
nj^a = ^c- ^, and U is the on-site Coulomb repulsion. 




FIG. 1: The t - t' Hubbard chain. 



The model can be viewed either as a single chain with 
both nearest- and next-nearcst-ncighbor hopping or, as 
illustrated in FigQ as a system of two coupled chains. 
The former view is appropriate for t ^ t' , the latter for 
t' > t. 

For t' = 0, we recover the ordinary Hubbard model 
which is exactly (Bethe Ansatz) solvablei^ In the case of 
a half-filled band, the ground state is insulating for arbi- 
trary positive values of U ; the charge excitation spectrum 
is gapped while the spin excitation spectrum is gaplessi^ 
For U <^ t the charge gap Ac is exponentially small, 
Ac w Vuie-'^''*/^ , while A^ w C/ for [/ > t& 

For t' ^ the model is no longer integrable except in 
the non-interacting limit, [/ = 0, where H is diagonal- 
ized by Fourier transformation and has a single-electron 
spectrum 

e(fc) = -2icosfc-|-2t'cos2fc. (2) 

For t' < 0.5t, the electron band has two Fermi points at 
kp ~ ±7r/2, separated from each other by the Umklapp 
vector q = TT (see FigEJ- In this case, a weak-coupling 
renormalization group analysisii predicts the same behav- 
ior as for t' — because the Umklapp term of order U is 
not modified; it again leads to the dynamical generation 
of a charge gap for U > 0, while the magnetic excitation 
spectrum remains gapless. 

In the strong-coupling limit, U ^ t,t', the charge sec- 
tor is gapped, while the spin sector can be mapped onto 
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FIG. 2: Single-particle dispersion relation of the t — t' chain 
for (a) t' = OAt and (b) t' = t. 



a frustrated Heisenberg chain 

H = J2iJS,- S,+i + J' S, • S,+2) (3) 
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with J = -it^/U and J' = W'^/U. This model has been 
extensively studied using a number of different analytical 
methodsSiiSSiS and has been found to develop a spin gap 
for J'/J ~ [f /tf > 0.24132i*2S and incommensurate 
antiferromagnetic order for J'/J > 0.5m^ This picture 
has been confirmed numerically^iii^ 

For t' > 0.5t, the Fermi level intersects the one-electron 
band at four points (±fc^) . This is the origin of more 
complex behavior for weak and intermediate values of U. 
For weak coupling [U <^t), the ground-state phase dia- 
gram is well understood in the two-chain limit [t' ^ t)& 
In this case, the Fermi vectors kp are sufficiently far from 
tt/2 to suppress first-order umklapp processes. Therefore, 
the system is metallic. The infrared behavior is gov- 
erned by the low-energy excitations in the vicinity of the 
four Fermi points, in full analogy with the two- leg Hub- 
bard modeli24 Thus, while the charge excitations are gap- 
less, the spin degrees of freedom are gappediSiiSiiiiSi^ 
Higher-order Umklapp processes become relevant for in- 
termediate values of U because the Fermi momenta fulfill 
the condition 'i{kp — kp) = 2n (at half-filling). There- 
fore, starting from a metallic region for small C/ at a 
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given value of t' {t' > 0.5t), one reaches a transition line 
U = Uc{t'), above which the system is insulating with 
both charge and spin gapsi^ The gross features of the 
phase diagram are depicted in Fig. |31 




FIG. 3: Qualitative sketch of the ground state phase diagram 
of the half-filled t — t' model. The solid line marks the metal- 
insulator transition. The dashed line indicates the transition 
from a gapless spin excitation spectrum at t' < 0.5t to the 
spin-gapped phase. 



III. BOSONIZATION 

A. Effective Theory 

We first consider the regime U,t' <^ t where bosoniza- 
tion is applicable. We linearize the spectrum in the vicin- 
ity of the two Fermi points kp = ±7r/2 and go to the 
continuum limit by substituting 



(4) 



where the operators ipRa{x) and ipLaix) are the right and 
left components of the Fermi field, respectively. These 
fields can be bosonized in a standard way^ 




where {4'Ra) are the right(left)-moving Bose fields 
and a is the infrared cutoff. We define 4>cr = 't'Ra + 't'La 
and introduce linear combinations, ipc — {(j)i+(j)i)V2 and 
(fis = {4>i — (j)i)V^, to describe the charge and spin degrees 
of freedom, respectively. Correspondingly, we introduce 



the conjugate fields 9c = 4>l<j — 4'Rit and = (0t+^i)v 2 
and i9s = (0| — 6i)^/2. After a simple rescaling, we arrive 
at the bosonized version of the Hamiltonian 



where both the spin part 



3(V87r(p5)|, 
and the charge part, 

n 

3(\/87r<^c)|, 



(5) 



■ cos 



(6) 



are described by the massive sine-Gordon model, with 
parameters 



{K^ - 1) -2m° = 2m° « -Vj-Kt . 



(7) 



There is an important difference between Hs and He due 
to the different stiffness constants. In the spin sector with 
Kg ~ 1, the system is in the weak-coupling limit and 
scales to a Gaussian model with gapless spin excitations. 
In the charge sector with Kc < 1, the system is in the 
strong coupling regime and the low-energy behavior is 
dominated by the cosine term. In the ground state, the 
field ifc is pinned at one of the minima of the cosine 
term and, correspondingly, there is a finite energy gap 
for charge excitations. 

Let us now discuss what happens when t' increases and 
reaches values of the order of t/2, where two additional 
Fermi points appear in the band structure. For spin 
degrees of freedom, new scattering channels appear at 
t' = t/2, and the system scales to strong coupling. There- 
fore, a spin gap is expected to open for t' > t/2, very 
much like in the case of two coupled Hubbard chains^ 

For the charge degrees of freedom, the situation is more 
complicated (and more interesting) because the charge 
gap blocks new scattering channels until t' is made suffi- 
ciently large so that additional states emerge beyond the 
gapped region. Thus, for t' slightly above t/2, the Fermi 
momentum changes without affecting the Umklapp pro- 
cesses. The Hamiltonian is still given by Eq. ©, but in 
order to allow for a change of particle number around 
the Fermi points, we have to add a topological (chemical 
potential) term 



where 



2t' 



V2^ 



2t' ^ 



dx dx^Pc 



for t' < 0.5t 
for t' > 0.5t . 



(8) 



(9) 



He + SHc is the standard Hamiltonian for the 
commensurate-incommensurate transition)2SiS& which 
has been intensively studied in the past using 
bosonizationi^ and the Bcthe ansatziSS 



B. Commensurate-incommensurate transition 

We now apply the theory of commensurate- 
incommensurate transitions to the metal-insulator 
transition as a function of t'. At fi^s = and K^. < 1, 
the ground state of the field (pc is pinned at 

(10) 



87r(y5c|0)o ~ 2?™ . 



The presence of the effective chemical potential makes it 
necessary to consider the ground state of the sine-Gordon 
model in sectors with nonzero topological charge. Using 
the standard expression for charge density in the case of 
two Fermi points^S^ 



Pc{x) 



1 



+ ^2fcF cos{2kFx) sin(\/27r(pc) cos(-\/27r(^s) 
-I- A^kp cos{Akpx) cos^VSnipc) , (11) 



we observe that the pinning of the field tpc in one of the 
minima 1)10(1 suppresses the 2kp charge fluctuations and 
stabilizes the Akp component. Any distortion of the ikp 
charge distribution would require an energy greater than 
the charge gap. This competition between the chemical 
potential term and conimcnsurability drives a continu- 
ous phase transition from a gapped (insulating) phase at 
Meff < Meff t° * gapless (metallic) phase at 



Meff > /"off = 



(12) 



where A^, is the charge gap at /ios = 0. 

We now separately consider the qualitative behavior of 
the system in the following three parameter regimes: (i) 
t' < OM, (a) OM < t' < t'^, and (Hi) t' > t'^, illustrated 
in Fig^ In regime (i), t' < 0.5t, we expect a charge 
gap Ac{U, t') w Ac(C/, t' = 0) and no spin gap, as in the 
simple Hubbard model {t' = 0). In regime (ii)^ Q.bt < 
t' < t'^, the spin gap opens while the charge gap is reduced 



A,((7,t') = A,(C/,0.5i) 



Mcff , 



(13) 



where /ieff is given by Eq. Therefore, the charge 

gap decreases with increasing t' and tends to zero at a 
qualitatively given by 



Ac(C/,0)-2t'^ + tV2ic = 0. 



(14) 



In regime (Hi), t' > t'^, the behavior of the system is char- 
acterized by four Fermi points, ±kp. The charge excita- 
tions are gapless, while the spin excitations are, generi- 
cally, gappediSi24 Charge fluctuations will be character- 
ized by two dominant periodic modulations with wave 
vectors 2k p and 2k p. For t' slightly larger than t'^, the 



usual charge-density wave {2kp ft 
a long-wavelength modulation at 



tt) is accompanied by 



£(k) 

(i) \ 










V < 0.5 







t' > 0.5 



M^eff > M^c 



2fc; 



(15) 



FIG. 4: Sketch of the energy dispersions for the three regimes 

(i), (a), and (Hi). 



As soon as the new set of states in the vicinity of the 
Fermi points ±fc^ appears, a channel for higher order 
umklapp scattering processes opens. In the sector of 
the phase diagram characterized by two-band behavior 
{f ^ t), these processes are responsible for the opening 
of a charge gap at the transition linei^ Therefore, in the 
parameter range where the renormalized one-band (Hub- 
bard) gap (|13|) becomes exponentially small, a crossover 
to the regime of two-band behavior takes place. There- 
fore, the linear decay of the charge gap as a function of 
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t' crosses over to exponential behavior. The evolution of 
the charge gap as a function of t' is sketched in Fig. [S] 





t'c 






FIG. 5: Sketch of the charge gap as a function of the param- 
eter t' . The inset shows an enlargement of the vicinity of the 
transition point. 



C. Two-chain limit 

To conclude our analysis, we discuss the limit of strong 
next-nearest-neighbor hopping {t' ^ t). For t — 0, the 
system is decoupled into two half-filled Hubbard chains 
and, for arbitrary [/ > 0, the ground state corresponds to 
a Mott insulator. The origin of the insulating behavior 
is the commensurability of umklapp scattering between 
the Fermi points, located at ±7r/4 and ±37r/4. When 
t this commensurability is lost. The Fermi points 
are shifted with respect to their values at i = 0, and the 
Fermi energy (the chemical potential for U ^ 0) moves 
away from to Ef « —t'^/2t' (for t <^ t'). For large 
enough values of t, the system is therefore expected to 
be metallic. 

In order to estimate the location of the Mott transition, 
we can use a similar argument to the one given above for 
t' K 0.5. As long as the chemical potential is smaller 
than the charge gap, the system remains an insulator. A 
transition to a metallic phase is expected to occur for e f 
of the order of A^, i.e., for t'^ w (Ut'^)^ cxp{-2nt' /U). 
A qualitative sketch of the phase diagram is given in Fig. 
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Mott 
Insulator 




t=0.5t' 



FIG. 6: Qualitative phase diagram of the half-filled repulsive 
t — t' Hubbard chain. A gapless charge excitation spectrum 
(metallic phase) exists ai U — for arbitrary t and t' and 
for ?7 > in the sector of parameter space below the "roof" 
covering the area U < Uc between the lines t' = 0.5t and t = 
in the U — plane. 



of lengths between L = 32 and L = 128 sites, keeping up 
to TO = 1000 density-matrix eigenstates. As we shall see 
in the following, the finite-size effects are quite large in 
certain parameter regimes, so that a careful finite-size 
scaling must be carried out. 



Transition line 




IV. NUMERICAL RESULTS 



In order to investigate the detailed behavior of the 
metal-insulator transition and to test the validity of the 
picture obtained from bosonization, we have carried out 
numerical calculations using the DMRGi^ We have cal- 
culated the properties of the ground state and low-lying 
excited states for systems with open boundary conditions 



FIG. 7: The metal-insulator transition line in the t — t' — U 
model with t = 1 obtained from DMRG studiesi^ (black 
circles) and from Eq. 11411 (solid line). 

The critical behavior of the metal-insulator transition 
as a function of U/t for t' > 0.5t can be obtained from 
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the behavior of the electric susceptibihty, which diverges 
in going from an insulator to a metalii^*2S In Fig. [7| we 
display the transition line in the t — t' — U model at t = 1 
obtained from the DMRGMiSa and from Eq. The 
agreement between the DMRG results and Eq. H14(l is 
remarkably good. 



B. Spin and charge gaps 

In order to investigate the predictions of the continuum 
theory, we calculate the charge gap, defined as 

1 



0.15 



A, = - [£;o(^ + 2,0) + Eo{N - 2, 0) - 2EoiN, 0)] 



and the spin gap, 



A, = Eo{N,l)~EoiN,0), 



(16) 



(17) 



where Eo{N, S) is the ground-state energy for N particles 
and spin S* on a chain of fixed length L, using the DMRG. 

In Fig. IHl we display the spin gap as a function of t' 
at U/t = 2 and U/t = 3 for various values of the chain 
length L. As can be clearly seen, for < t' < 0.5t the 
spin excitation spectrum on finite chains does not depend 
on t' . For t' < 0.5t, the value of the spin gap is found 
to coincide with that of the half-filled Hubbard model 
(t' — 0) which vanishes in the infinite-chain limit [see 
Fig.|i(a)]. 

A clear change in the t'-dependence of the spin gap at 
U/t = 2 takes place at t' = O.bt, indicating the devel- 
opment of a new phase in the spin sector. It is known 
from other studies^tiSiiiiSiSi that a spin gap opens at a 
critical value of t' which is is approximately at or slightly 
above t' = 0.5i, becoming weakly larger at intermediate 
U values. 

In Fig. Eta) , we display the spin gap plotted as a func- 
tion of the inverse chain length for three values of t' near 
the transition at U/t = 3. At t' = 0.55t, the spin gap 
clearly scales to zero at infinite system size, with the val- 
ues at a particular system size virtually identical to the 
t' = case and the scaling predominantly linear in inverse 
system size. For t' = 0.6t and t' = 0.65t, the dominant 
scaling term is quadratic rather than linear in 1/L and 
there is clearly scaling to a finite value of the gap. For 
t' ~ 0.65t, the size of the extrapolated gap is smaller than 
for t' = 0.6t, and there is a slight upturn in the gap at the 
largest system size, which, however, is not significantly 
larger than the estimated error of the DMRG calcula- 
tion, approximately the symbol size. However, for larger 
values of t', the finite-size behavior becomes less regular, 
as can be seen in Figs. Ela) and (b). This behavior is 
due to the appearance of an incommensurate wave vec- 
tor characterizing the spin excitations that occurs when 
a substantial density of states at all four Fermi points 
develops and makes it virtually impossible to carry out 
a well-controlled finite-size scaling for larger values of t'. 

The transition associated with the opening of the spin 
gap is independent from the insulator-metal transition, 
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FIG. 8: Spin gap of as a function of t' for (a) U/t = 2 and 
(b) U/t = 3. 



as can be clearly seen for U/t = 3 (where the effect of 
fluctuations is reduced). As is shown in Fig. ^h) and 
Fig. inia), the spin gap opens for t'^ > 0.55t, while the 
insulator-metal transition takes place at t'^ ~ 0.65i (see 
below) . Note that the critical value of the next-nearest- 
neighbor hopping amplitude, corresponding to an open- 
ing of the spin gap at U/t = 3, t'^ > 0.55i deviates from 

the line t'^ > O.bt. Our findings agree with previous 
studiesi«-i^i",i5,24 

A plot of the finite-size extrapolation of the charge gap 
is displayed in Fig.lHJb) for various values of t' . As can 
be seen, the behavior is well-behaved for values of t' from 
to 0.8. For < t' < 0.6t, the scaling has a substantial 
positive quadratic term in 1/L and the gap is finite. For 
t' — 0.65t and 0.8t, the extrapolated gap clearly vanishes 
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and there is a negligible or negative quadratic contribu- 
tion. For t' > 0.8, the finite-size effects become irregular 
due to incommensurability of the charge excitations, and 
finite-size extrapolation becomes difficult. 
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FIG. 9: (a) Spin gap and (b) charge gap of as a function of 
I/L for U jt — 3 and various values of t' /t. 

In Fig. EH the L = oo extrapolated value of the charge 
gap is displayed as a function of t' for U/t = 2 and U/t = 
3. There is a clearly defined insulator- metal transition at 

= 0.55t at U/t = 2 and tc = 0.65t for U/t = 3. Note 
that the charge gap goes smoothly to zero above t' = 
0.5t for U/t = 3. The inset in Fig. If 01 shows the charge 
gap for J7/t = 3 as a function of the parameter ^^s = 
2t' - t^/2t' for 0.5i < t' < 0.85t As can be seen, the 
charge gap drops off approximately linearly with ^eff, in 
agreement with Eq.^l For U/t = 2, there is a somewhat 
irregular behavior of the charge gap near the t' = 0.5t. 



In particular, there is a small peak exactly at t' = 0.5. 
The finite-size scaling for this point is completely regular, 
however, and we estimate the size of the total error, due 
to both the extrapolation and the DMRG accuracy, to be 
less than the symbol size. Therefore, in our estimation, 
the peak at t' = 0.5 is a real effect. For <' = 0.55, the 
value of the extrapolated charge gap is slightly below 
zero. This is due to errors in the finite-size extrapolation 
due to slightly irregular behavior with system size. 



j_i i_i i_j i_j I ^ I \ I ^ I 




0.2 0.4 0.6 0.8 

t'/t 

FIG. 10: Charge gap of as a function of t' for U/t = 3 (black 
circles) and U/t = 2 (open circles) . The inset shows the 
charge gap as a function of the parameter /Xeff for 0.5 < t' < 
0.85t. 



C. Spin and charge densities 

Valuable insight into the nature of the insulator-metal 
transition can be obtained by studying the charge and 
spin density distributions. 

We first examine the local charge density {n/) in the 
ground state with L = 64 sites and U/t = 3 as i' is varied. 
At t' = 0.6t, as can be seen in Fig. Illf a). the commen- 
surate charge distribution characterizing the insulating 
phase (i.e., (n^) = 1) is reached within a few lattice sites 
from the edge. The boundary effect is relatively weak and 
short-range. The insulator-metal transition manifests it- 
self at t'^ = 0.65t via the appearance of incommensurate 
modulations in the charge distribution. One can clearly 
see the long-wavelength modulations in density originat- 
ing from the opening of the Fermi surface for |fc| < Kp. 
At t' = O.St, the incommensurability increases and the 
presence of long wavelength modulations of the charge 
density become more evident. 

In Fig. llir b). we plot the charge distribution at t' — 0.8 
for various values U . In the metallic phase, at U/t — 2 
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FIG. 11: Charge distribution {rn) in the ground state of 
the t - t' - U chain with (a) L = 64, U/t = 3 and t' = 
0.6t, 0.65t, 0.8t and (b) for t' = 0.8t and U/t = 2, 3, 6, 8. 



and U/t ~ 3, the charge distribution pattern is strongly 
incommensurate. At J7/t = 6 and U/t = 8, the well- 
established commensurate pattern of the charge distri- 
bution characterizes the insulating phase. 

Let us now consider the spin density distribution. As 
has already been stressed above, due to spin-charge sep- 
aration, the magnetic degrees of freedom in the t — t' — U 
Hubbard model are not influenced by the destruction 
of the "holon Fermi surface" caused by the dynamical 
generation of a charge gap. Instead, the development 
of the incommensurate spin distribution in the insulat- 
ing phase is completely determined by the competition 
between the nearest-neighbor and next-nearest-neighbor 
spin exchange interactions and reflects the itinerant na- 



FIG. 12: Spin distribution l^S^ in the Sjotai = 1 state of 
the frustrated Heisenberg chain with L = 64 and J' / J = 
0,0.4,0.5, and 0.6. 



turc of the model only when the system is very close to 
the insulator-metal transition or is in the metallic phase. 
To see this, we examine the behavior of the spin distri- 
bution {Sf) in the S^^^ = 1 excited state of the t — t' — U 
Hubbard model, comparing it to the effective model for 
strong coupling, the frustrated Heisenberg chain. 

We first examine the limiting case of the frustrated 
Heisenberg chain, Hamiltonian In Fig. ^1 we plot 
(5f) calculated using the DMRG in the Sf^^ = 1 state 
for various values of the next-nearest-neighbor exchange 
J' . At J' = 0, the unfrustrated case, there is one node 
in the distribution of the spin density, which corresponds 
to two spatially separated spin-S* =1/2 spinonsi^ These 
two spinons characterize the spin excitations up to the 
Majumdar-Ghosh point J' = 0.5 J {f ~ 0.7t)^ For 
larger values of the next-nearest-neighbor exchange, the 
spin distribution becomes incommensurate In the plot 
for J'/ J — 0.6, the absence of the commensurate two- 
spinon structure of the excitations and the appearance 
of the incommensurate wave vector can be seen. 

We next examine the behavior for the t — t' — U chain. 
In Fig ^1 we show the spin distribution (5f ) in the 
•^totai ^ 1 state with L — 64, U/t — 3 and various val- 
ues of t' . The single- node pattern characteristic of two 
spinons, as found in the unfrustrated and weakly frus- 
trated Heisenberg chain, can be seen for t' = 0.6t and 
smaller. At t' = 0.7t an incommensurate pattern ap- 
pears. This is just above the value of t' at which the 
charge gap goes to zero, t'^ ~ 0.65i, (see Fig. and 
the system becomes metallic. However, for t' > 0.7t the 
point at which the spin-density distribution becomes in- 
commensurate is independent of the insulator-to-metal 
transition. To demonstrate this, we plot the spin density 
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FIG. 13: Spin distribution (Sf ) in the ^totai = 1 state of the 
t-t' -U chain with L = 64, U/t = 3 and t' = 0, OAt, 0.6t, 
and 0.7t. 



distribution in the S^^^^i = 1 state for increasing on-site 
repulsion U for t' = 0.6t and t' = 0.8t in Fig. [H As 
can be seen in Fig. I14f a'). for t' ~ 0.6t the spin-density 
distribution is incommensurate in the metallic phase at 
U/t = 2, but acquires the commensurate single-node pat- 
tern characteristic of two spinous for on-site repulsions 
oi U/t = 3 and higher. In contrast, for t' = 0.8t, Fig. 
I14f b'l. the incommensurate pattern of the spin distribu- 
tion in the metallic phase (U/t = 2) remains not only in 
the vicinity of the insulator-metal transition at U/t = 3 
but also deep into the insulating phase at U/t = 8 and 
U/t = 10. 
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D. Excitations for Large t' 

We examine the behavior of the t—t' — U chain for large 
next-nearest hopping [t' ^ i) , a limit which corresponds 
to two chains coupled with a weak zigzag hopping. In 
particular, we numerically investigate the transition from 
a two-chain (four-Fermi-point) metallic regime at weak 
U to the strong-coupling regime, for which the effective 
model is two spin-5 = 1/2 Heisenberg chains coupled 
with a frustrating zigzag interaction at [/ 3> i' 3> i, i.e., 
J' > J. 

In Fig. 1151 we show the charge distribution {rii) in the 
ground state of the t—t' — U chain for t' = 3t, L = 64, and 
U/t = 4,8, 10, and 20. In the metallic phase {U/t = 4 
and U/t — 8), the relative charge density on the two 
chains alternates with relatively strong amplitude and 
is modulated by an incommensurate wave vector. As 
U/t is increased to 10, the amplitude of the charge fluc- 
tuations between the chains is strongly suppressed and 



FIG. 14: Spin distribution (Sf) in the Stotai = 1 state of 
the t — t' — U chain with L = 64 and for (a) t' = 0.6t with 
U/t = 2, 3, 6, 8 and (b) t' = 0.8t with U/t = 2, 4, 8, 10. 

the incommensurate structure of the charge distribution 
changes to a very distinctive two-bubble pattern. This 
two-bubble distribution corresponds to two spatially sep- 
arated electron-hole pairs with an electron on even chain 
and a hole on odd chain in the left part of the lattice 
and the inverse configuration in the right half of the lat- 
tice. (Note that the ends break the symmetry between 
the two chains.) Deeper into the insulating phase, at 
U/t — 20, the charge density is almost smooth and equal 
between the chains, except for a small residual alterna- 
tion of charge density due to the end effects. 

In Fig. El we show the spin density distribution (5f ) 
in the S^^^^^ = 1 state for t' = 3t and U/t = 4, 8, 10, 20. 
As can be seen, the incommensurate pattern of the 
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and spin gap with increasing next-nearest-neighbor hop- 
ping ampHtude t' . We have found evidence for a spin 
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FIG. 15: Charge distribution (rii) in the ground state of the 
t-t' ~U clmin with L = 64, t' = 3t and U/t = 4, 8, 10, 20. 
Black circles correspond to odd and open circles to even lattice 
sites. 



spin distribution in tlic metalhc phase at U/t — 4 and 
U/t = 5 transforms in the insulating phase {U/t ~ 10 
and U/t ~ 12) into a pattern which, for each chain, is 
similar to that of the single Heisenberg chain in the com- 
mensurate phase. This pattern, seen in, e.g.. Fig. 1121 for 
J'/ J ~ 0, 0.4, 0.5, corresponds to two-spinon excitations, 
as previously discussed. This indicates that the system 
behaves as two identical, weakly coupled 5* = 1/2 Heisen- 
berg chains at large U. 



V. CONCLUSIONS 

We have carried out a combined analytical and numer- 
ical analysis of the insulator-metal transition in the half- 
filled one-dimensional t~t'~U model. Using the weak- 
coupling bosonization approach, we have shown that the 
gross features of the transition from an insulator to metal 
as a function of next-nearest-neighbor hopping t' can be 
described within the standard theory of commensurate- 
incommensurate transitions. We have derived an explicit 
expression for the critical line t'^{U) separating the metal- 
lic phase from the spin-gapped insulator. We have also 
argued that in the vicinity of transition line a crossover 
to infinite-order critical behavior occurs. 

Using DMRG calculations on chains of up to L = 128 
sites, we have performed a detailed numerical analysis of 
the excitation spectrum and the charge and spin density 
distributions in various sectors of the phase diagram. In 
particular, we have studied the evolution of the charge 



FIG. 16: Spin distribution (Sf ) in the Stotai = 1 state of the 
t-t' -U chain with L = 64, t' = Zt and U/t = 4, 8, 10, 20. 
Black circles correspond to lattice odd and open circles to 
even lattice sites. 



gap in the parameter range 0.5i < t' < t'^, in agreement 
with previous studies. We have shown that the change 
in the topology of Fermi surface at the insulator-metal 
transition is reflected in the appearance of incommensu- 
rate modulations of the charge density. Incommensurate 
spin-density distributions in the triplet sector are always 
present in the metallic phase, but can also appear inde- 
pendently in the spin-gapped insulator due to frustration. 

For t' ^ t, we have argued that the insulator-metal 
transition can be best understood starting from the limit 
of two uncoupled chains. At small U, turning on the 
zigzag coupling between the chains destroys the com- 
mcnsurability present for a single chain, and leads to 
a metallic phase. At large U, the system is insulating 
and behaves as two weakly coupled Heisenberg chains. 
We have estimated that the insulator-metal transition in 
this regime occurs when the shift in the Fermi energy is 
comparable to the size of the charge gap in the isolated 
Hubbard chain. 
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